■ 


Pra-lact  No.  53 
Task  3  Technical 


SIMULATION  OF _ADAPTIVE _£LOSED-LOOP  J10NTRGL__PROCEDURES 
FOR  ELECTRONIC  COUNTERMEASURES  . 


Sponsored  by: 

DEFENSE  ADVANCED  RESEARCH  PROJECTS  AGENCY 


ARPA  Order  No.  3613 
Monitored  by  NAVELEX  under 
Period  of  Contract:  27  Decei 
Reporting  Period:  27  Decei 


RGia 


SjARfb  Sr  Jr 


Approved  by :  . 

Q.  f.  &  ■ ' 

J.P.  Marsh 
Manager , 

Control  Technology  Program 


Prepared  by: 

A.  Joseph i  Rccteore 
Paul  L./' Cowell 
Scott/Foster 


:  DISTRIBUTION  STATEMQfT  A  ; 

Approved  for  public  release; 

Di-.lithudon  Unlimited 

The  views  and  conclusions  contained  in  this  document  are  those  of  the 
authors  and  should  not  be  interpreted  as  representing  the  official 
policies,  either  expressed  or  implied,  of  the  Defense  Advanced  Research 
Projects  Agency  or  the  U.S.  Government." 

yir1333  J4»80  9  26  029 


SMTIMS  CONTROL,  NC  ■  1801  PAGE  MILL  ROAD  □  PALO  ALTO,  CA  84304  ■  TELEX  348-433  ■  <4I5>  494-1165 


L 


TABLE  OF  CONTENTS 


1.  INTRODUCTION .  1-1  | 

1.1  Description  of  the  Simulation .  1-2  1 

1.2  Coordinate  Systems  .  1-6  * 

1.2.1  Topocentric-Horizon  (SEZ)  Coordinates  .  1-6 

1.2.2  Antenna  Coordinate  System  and  Beam  Geometry  .  .  .  1-6 

1.2.3  Target  (RPY)  Coordinate  System  and 

Target  Geometry  .  1-18 

1.3  Radar  Signal  Processing  .  1-20 

1.3.1  Digital  Filters .  1-21 

1.3.2  Automatic  Gain  Control .  1-23 

1.3.3  Antenna  Servomechanism  Models  .  1-24 

1.3.4  Error  Detector  in  Conical  Scan  System  .  1-25 

1.4  The  Scenario .  1-26 

1.4.1  Target  Model .  1-26 

1.4.2  Target  Trajectory  .  1-31 

2.  CONICAL  SCAN  SYSTEMS .  2-1 

2.1  Conical  Scan  Threat  Tracking  Loop .  2-1 

2.2  Simulation  of  ECM  Techniques .  2-5 

2.2.1  Parametric  Controller  Against  COSRO  .  2-7 

2.2.2  Input/Output  Controller  Interface  for  COSRO  .  .  .  2-12 

3.  PS UEDO- MONOPULSE  SYSTEMS  .  3-1 

3.1  Psuedo-Mbnpulse  Tracking  Loop  .  3-1 

3.2  Simulation  of  ECM  Techniques .  3-4 

3.2.1  Parametric  Controller  Against  Psuedo-Monopulse  .  .  3-4 

3.2.2  Input /Output  Controller  Interface  .  3-6 

4.  MONOPULSE  SYSTEMS  .  4-1 

4.1  Monopulse  Threat  Tracking  Loop .  4-1 

4.2  Simulation  of  ECM  Techniques .  4-5 


I 


» 


TABLE  OF  CONTENTS  (Cont’d) 


5.  SAMPLE  RUNS .  5-1 

I  5.1  Sample  Runs  Against  COSRO .  5-2 

5.2  Sample  Runs  for  Psuedo-Monopulse .  5-25 

5.3  Sample  Runs  for  Monopulse .  5-32 

REFERENCES .  R-l 

> 


t 


ill 


A 


1 


INTRODUCTION 


^Techniques  to  accomplish  the  goals  of  electronic  countermeasures  (ECM) 
have  been  developed  in  the  past  by  exploiting  a  specific  feature  of  a  sys¬ 
tem  in  a  heuristic  manner.  In  general,  emphasis  has  been  placed  upon  the 
development  of  a  technique  which  has  some  effectiveness  against  a  particu¬ 
lar  threat  system,  and  the  choosing  of  the  parameters  of  the  technique  has 
received  little  attention.  This  report  documents  the  results  of  the  third 
task  of  the  work  Systems  Control,  Inc.  (SCI)  has  performed  to  apply  advanced 
concepts  from  modern  estimation  and  control  theories  to  provide  improved 
performance  of  ECM  systems.  In  this  program,  approaches  to  the  design  of 
techniques  for  ECM  are  developed,  especially  as  related  to  the  control,  or 
time-varying  choice  of  parameters  of  the  technique,  so  as  to  achieve  a 
specified  performance  objective.  <5- - 

Task  I  Technical  Report  [1]  is  concerned  with  the  development  of 
control  procedures  to  vary  the  parameters  of  the  ECM  technique.  These 
procedures  are  derived  from  the  disciplines  of  control  and  estimation 
theories.  They  fall  into  the  category  of  adaptive  closed-loop  jamming, 
where  observables  related  to  the  radar  (e.g.,  the  time-varying  power  of  a 
conical  scan  system,  the  actual  threat  antenna  pointing  angle,  etc.)  are 
used  in  a  feedback  structure  to  vary  the  control  variables  of  an  ECM 
technique  (e.g.,  the  amplitude  of  the  pulses  in  a  scan  rate  modulation 
ECM,  the  relative  amplitude  and  phase  of  a  cross-eye  ECM,  etc.).  The 
unknown  nature  of  the  threat  radar  is  what  calls  for  adaptive  controllers. 

Task  II  Technical  Report  [2]  is  concerned  with  the  application  of  the 
techniques  described  in  [1]  to  three  particular  threat  systems:  conical 
scan  on  receive  only  (COSRO),  pseudo-monopulse,  and  monopulse.  Each 
threat  system  is  analyzed  in  terms  of  available  observables  and  an  appropriate 
ECM  technique.  The  connection  between  these  two  is  provided  by  the  adaptive 
controller,  which  takes  the  observable  and  feeds  it  back  to  the  ECM  control 
variable.  The  form  of  the  controller  takes  on  one  of  two  forms,  described 
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in  [1].  The  first  is  a  parametric  controller,  which  exploits  the  actual 
structure  of  the  threat  radar  and  controls  using  ideas  from  stochastic 
approximation.  The  second  is  an  input/output  controller,  which  assumes  an 
auto- regressive  moving  average  (ARMA)  model  relating  the  observables  to 
the  control  variables,  and  does  a  model  fit  and  control  based  on  the  model. 

Both  of  these  approaches  are  employed  against  the  specific  systems  of  interest. 

The  current  report  is  concerned  with  the  simulation  of  the  threat/ 
target  interaction  to  verify  and  explore  aspects  of  the  performance  of  the 
adaptive  closed-loop  approach  to  ECM.  A  scenario  is  postulated  in  which 
the  threat,  which  is  a  ground-based  or  missile  radar,  illuminates  a  tar¬ 
get,  which  is  an  aircraft  or  a  ship.  The  emphasis  in  the  current  report 
has  been  on  the  problem  of  defending  an  aircraft  against  a  ground-based 
threat,  but  elements  of  the  antiship  missile  defense  problem  are  included. 

The  latter  problem  will  receive  more  emphasis  in  the  following  phases  of 
the  project.  Models  used  and  results  obtained  are  described. 

This  report  assumes  a  faiily  high  degree  of  familiarity  with  the 
previous  task  technical  reports.  Thus  details  of  the  ECM  techniques  and 
derivation  of  the  control  procedures  are  omitted.  Detail  related  specifi¬ 
cally  to  the  computer  simulation  are  included,  as  well  as  results. 

1.1  DESCRIPTION  OF  THE  SIMULATION 

The  purpose  of  the  simulation  is  to  compare  and  evaluate  the  perfor¬ 
mance  of  various  ECM  systems  designed  to  defeat  one  or  more  types  of 
tracking  radars.  The  simulated  system  is  composed  of  two  subsystems,  the 
threat  and  the  target.  The  threat  system  is  a  tracking  radar  which  actively 
illuminates  a  target  to  establish  its  position.  The  target  is  the  aircraft 
or  ship  being  tracked  and  its  ECM  and  ESM  equipment. 

A  system,  which  may  be  defined  as  a  collection  of  entities  in  a  regu¬ 
lar  interaction,  is  modeled  by  a  body  of  selected  information  which  describe 
the  salient  system  features  to  be  studied.  In  general,  the  models  are 
specified  by  the  four  structural  elements  described  below. 


The  entities:  objects  which  comprise  the  system. 

The  system  is  composed  of  smaller  subsystems  and  the  correspond¬ 
ing  entities  may  be  grouped  into  a  hierarchy  of  more  general 
entities,  such  as: 

Threat 

Airframe  (if  missile  or  aircraft)  dynamics 
Radar 

Antenna 
Servo  motors 
Tracking  loop 
Detector 

Signal  processing 
Servo  motor  driver 

Target 

Airframe  (if  aircraft)  or  ship  dynamics 
E W  equipment 

ESM  antenna 
Sensor  processing 
ECM  controller 
ECM  technique  generator 
ECM  antenna. 

Each  of  the  threat  and  target  systems  is  composed  of  four  basic 
subsystems : 

•  a  sensor  to  detect  radar  energy 

•  a  signal  processor  and  controller 

e  a  transmitter  to  generate  radar  energy 
e  a  reflector  of  incident  radar  energy. 

The  attributes:  properties  of  the  entities. 

The  attributes  describe  various  properties  of  the  entities, 
such  as  locations,  pointing  directions,  antenna  and  filter 
parameters,  and  radar  cross  section. 

The  activities:  causes  of  change  in  the  system. 

The  target  and  threat  may  be  thought  of  as  two  separate  but 
interacting  systems  which  have  a  spatial,  a  temporal,  and  a 
physical  relation.  The  two  systems  are  spatially  related  in  a 
coordinate  system.  Each  may  move  so  that  the  relation  may  be 
altered.  The  temporal  relation  is  constrained  to  synchronous, 
pulse-by-pulse  operation.  The  simulation  is  event-oriented  in 
that  it  follows  a  history  of  activities  which  are  applied  to 
different  entities. 

This  representation  of  the  target-threat  interaction  leads  to 
restrictions  on  the  types  of  relationships  that  can  be  modeled. 

The  timing  in  the  current  simulation  is  based  on  an  interval- 
oriented  clock  that  runs  at  the  pulse  repetition  frequency  (p.r.f.)  of 
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the  tracking  radar.  Because  events  occur  only  at  the  time  of 
the  simulation  clock  update,  the  threat  and  target  run  synchron¬ 
ously.  This  is  not  a  necessary  restriction  on  the  simulation, 
but  asynchronous  operation  has  not  yet  been  required  in  evaluat¬ 
ing  the  ECM  systems  developed  to  date. 

The  physical  relation  is  the  transmission,  reflection,  and 
reception  of  radar  energy  by  both  systems.  The  basic  activities 
which  alter  the  entities  are: 

•  Transmit  a  tracking  pulse 

•  Transmit  an  ECM  pulse 

•  Generate  a  skin  return 

e  Detect  the  ECM  and  skin  return 

•  Move  the  airframe  or  ship. 

Some  of  these  activities  are  deterministic,  such  as  pulse  trans¬ 
mission,  while  others  are  stochastic;  the  outputs  are  not  completely 
specified  in  terms  of  the  inputs.  Clearly,  these  activities  may 
be  broken  down  even  further  into  more  fundamental  activities  which 
change  the  most  basic  entities  in  the  system. 

(4)  The  environment:  changes  which  affect  the  system  but  occur 
outside  of  it. 

The  system-environment  boundary  includes  the  whole  simulation; 
in  other  words,  there  are  no  exogenous  activities,  those  that 
lie  outside  the  system. 

The  simulation  is  structured  as  shown  in  Figure  1-1  to  model  a  single 
engagement  between  a  target  airframe  or  ship  moving  on  a  pre-specified 
trajectory  past  a  stationary,  ground-based  or  moving  threat.  The  threat 
motion  can  be  in  response  to  the  radar  return  and  ECM  of  the  target,  such 
as  will  be  the  case  in  a  missile  threat.  An  engagement  is  divided  by  the 
user  into  a  number  of  segments  by  specifying  the  length  of  time  which 
will  elapse  in  each  segment,  where  time  is  measured  in  pulses.  Inputs 
to  the  simulation  define  the  target  and  threat  parameters,  the  target  and 
threat  trajectories,  and  control  the  simulation  itself.  Outputs  include  a 
complete  list  of  all  input  parameters,  aid  to  accurate  numerical  analysis 
or  debugging,  and  graphical  output  showing  temporal  sequences  or  histograms. 
These  outputs  are  under  the  control  of  a  user  who  determines  the  number  of 
pulses  occurring  at  the  threat  radar  p.r.f.  in  the  next  segment  of  the 
engagement  and  what  output  will  be  generated  for  it. 


The  major  part  of  the  modeling  of  the  threat,  target,  and  ECM  occurs 
in  the  tracking  loop,  shown  in  Figure  1-2.  (The  acquisition  or  search 
mode  is  essentially  the  same  as  the  tracking  mode  except  that  no  ECM  is 
allowed  and  the  antenna  is  forced  to  point  near  the  target  to  allow  the 
tracking  loop  to  stabilize.)  In  the  tracking  loop,  a  pulse  transmission 
from  the  threat  radar  illuminates  the  target.  The  resulting  skin  return, 
along  with  any  ECM,  are  detected  by  the  threat  to  provide  an  error  signal 
which  points  the  antenna  through  its  servo  motors.  In  order  to  accomodate 
a  variety  of  threat  antenna  and  signal  processing  components,  a  variety 
of  target  models,  and  a  variety  of  ECM  techniques,  the  tracking  loop  is 
broken  into  a  set  of  interacting  modules,  as  shown  in  Figure  1-3. 

1.2  COORDINATE  SYSTEMS 


The  physical  relationships  between  the  target  and  threat  as  they  move 
with  time  must  be  referenced  to  a  particular  coordinate  system.  This 
coordinate  system,  described  in  Section  1.2.1,  is  a  topocentric-horizon 
coordinate  system  called  the  SEZ  system,  for  South-jJast-2:  axis.  It  is 
thus  a  Cartesian  coordinate  system. 

The  antenna  gain  patterns  for  all  antennas  must  be  referenced  to  a 
translated  Cartesian  coordinate  system,  denoted  the  S'E'Z'  system.  This 
system  is  described  in  Section  1.2.2. 

Since  the  target  may  be  an  extended  target  (as  in  the  ship  case)  or 
its  ECM  technique  may  depend  on  spatially  separated  emitters  (as  in  the 
cross-eye  case) ,  a  separate  coordinate  system  must  be  defined  for  the 
target.  This  system,  described  in  Section  1.2.3,  is  called  the  RPY  system, 
for  Roll  axis-Pitch  axis-Yaw  axis.  It  is  thus  also  a  Cartesian  coordinate 
system. 

1.2.1  Topocentric-Horlzon  (SEZ)  Coordinates 

The  threat  and  the  target  are  located  in  a  common  coordinate  system 
fixed  on  the  surface  of  the  earth  called  the  topocentric-horizon  system, 
shown  in  Figure  1-4. 


Figure  1-3.  Modular  Structure  of  Simulation 


1-8 


Figure  1-4.  Topocentric-Horizon  Coordinate  System 


Figure  1-5.  Angles  Shown  for  Antenna  Pointing  Due  North 


-i 


The  origin  is  located  at  the  position  of  the  fixed,  ground-based  radar 
associated  with  the  threat,  or  at  an  arbitrary  place  if  there  is  no  fixed 
threat  radar.  The  S-E  plane  coincides  with  the  horizon  plane,  with  S 
pointing  south  and  E  pointing  east.  The  Z-axis  points  toward  the  Zenith. 
Any  vector  p  may  be  expressed  in  these  coordinates  by 

P  -  PgS  +  P£E  +  P£Z  . 

The  azimuth  and  elevation  angles  and  the  magnitude  of  the  vector  are  defined 
by 


AZ  =  arctan  pE/(-pg) 


E l  =  arcsin  p£/p 


o  -  <p*  ♦  4*  o\) 1/2 


The  inverse  transformation  is  given  by: 


Ps  *  -p  cos  E£  cos  AZ 
p£  =  p  cos  El  sin  AZ 
Pz  =»  p  sin  El  . 


1.2.2  Antenna  Coordinate  System  and  Beam  Geometry 

In  general,  an  antenna  gain  pattern  associated  with  the  antenna  will 
not  necessarily  be  aligned  with  the  antenna  axis.  A  coordinate  system 
attached  to  the  antenna  to  describe  the  orientation  of  its  gain  pattern 
is  shown  in  Figure  1-5.  It  has  translated  coordinates  S'E'Z',  correspond¬ 
ing  to  translated  South,  translated  East,  and  translated  Zenith.  The 
coordinate  system  axes  coincide  with  the  topocentrlc-horizon  system  axes 
when  the  antenna  is  pointed  due  north  (zero  azimuth  and  zero  elevation) . 
The  gain  pattern  location  is  described  by  its  squint  angle,  $,  and 
phase  angle,  y. 
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The  SEZ  coordinate  system  may  be  transformed  by  two  successive  rota¬ 
tions  into  the  S'E'Z'  coordinate  system,  as  shown  in  Figure  1-6.  The 

A 

first  rotation  is  made  about  the  Z-axis  through  angle  AZ,  the  azimuth 
angle.  The  second  is  made  about  the  new  E-axis  through  EJt,  the  eleva¬ 
tion  angle.  Denote  an  arbitrary  vector  by 

P  »  PSS  +  PeE  +  PZZ  -  p^S'  +  PgE’  +  p^Z’  . 

Then  the  coordinate  transformation  may  be  written 


When  the  antenna  axis  is  pointed  at  (Ao,Eq),  the  vector  along 

the  gain  pattern  axis  can  be  expressed  in  the  topocentric-horizon  system 
by 


r 

c 


R  (+A  )  R  (+E  )  r'  . 
z  o  y  o  c 


The  simulation  is  capable  of  modeling  antennas  with  multiple  beams 
whose  locations  are  defined  relative  to  the  bores ight  of  the  antenna  in 
the  antenna  coordinate  system. 

Figure  1-7  shows  the  parameters  used  to  define  the  beams,  each  of 
which  is  squinted  relative  to  the  antenna  boresight  by  an  angle  8^  and 
is  located  at  a  true  phase  angle  y^  from  the  elevation  axis.  Table  1-1 
indicates  the  names  of  the  quantities  used  in  the  simulation. 
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Table  1-1.  Beam  Geometry  Description  Input  to  Program  IANTG 


SYMBOL 

VARIABLE 

NAME 

DESCRIPTION 

N 

NB 

Number  of  beams 

8J 

BETA 

Squint  angle  of  beam 

vi 

GAMMA 

Initial  phase  angle  of  beam 

Gj  (0) 

GNAXIS 

Beam  gain  (if  transmitter) 

9j 

o 

BWDEG 

3  dB  (half)  beamwidth  in  degrees 

RCVARE 

2 

Beam  effective  aperture  in  m  (if 
receiver) 

When  the  phase  angle  of  any  beam,  Y^  ,  depends  on  time,  the  antenna 
system  is  described  as  a  conical  scanning  system. 

Antennas  with  conical  scanning  beams  are  modeled  by  computing  new 
beam  positions  in  the  antenna  coordinate  system  on  every  pulse. 

The  inputs  required  for  conical  scanning  antenna  system  are  shown 
in  Table  1-2.  For  each  pulse,  the  phase  of  each  beam  is  updated  in  the 
program  CSC,  if  the  beam  is  rotated  by  a  conical  scan  mechanism  according 
to 


Yk  -  (Y^j+dy)  mod(2rr)  . 

A  unit  vector,  r^,  along  each  beam  axis  in  the  antenna  coordinate  system 
is  computed  from 

r^  *  (-cos  6^,  -sin  3^  sin  y^,  +sin  3^  cos  Yk)T  • 
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Table  1-2.  Conical  Scanning  Antenna  Parameters  in 

Program  ICSG 


SYMBOL 

VARIABLE 

NAME 

DESCRIPTION 

f 

r 

PRFREQ 

Pulse  repetition  frequency 

f 

c 

CSFREQ 

Conical  scan  frequency 

*c 

CSPHAS 

Reference  signal  phase  shift 

dY 

DELGAM 

Conical  scan  phase  increment  per  pulse 

-  2tt  f  f_1 
c  r 

All  of  the  beams  rotate  together  in  a  conical  scan  system  about  the  antenna 
boresight,  maintaining  their  relative  phase  angles  constant.  The  reference 
signals  are  computed  relative  to  the  position  of  the  first  beam  by 

see*  f 

\  -  -sin(yk  +  *c) 

Ek*f  -  +cos(Yk  +  <1>C)  • 

In  the  tracking  radar,  the  servo  motors  drive  the  antenna  in  a  direc¬ 
tion  based  on  the  error  signal  derived  from  the  radar  returns.  Those 
corrections,  dA^  and  dE^  obtained  from  the  servo  motor  models,  are 
applied  to  the  current  antenna  position  and  the  true  position  of  the 
target  within  the  antenna  beams  is  obtained. 

The  servo  output  updates  the  antenna  pointing  angles: 

*k  •  Vi  +  %  ftl 

\  -  Vi +  dEk  *;1  • 
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The  current  antenna  boresight  axis  is  then  defined  by 


T 

b^  *  (-cos  cos  A^,  +cos  sin  A^,  sin  E^) 

Denoting  by  pk  the  vector  from  the  threat  to  the  true  target  position, 
the  true  azimuth  and  elevation  of  the  target  are  given  by 

tan  a£  -  -Pk  •  E(pk  •  S)  1 

sin  e£  -  pfc  •  Z  (pj-1  . 

The  beam  axes  are  rotated  into  the  SEZ  coordinate  system  by  the  suc¬ 
cessive  rotations 

sin  E°\ 

0  K 

cos  Ek  / 

The  angle  of  the  line-of-sight  to  the  target  from  each  beam  axis  is  then 
given  by 


coa  -  pk  •  rk 


,-l 


These  calculations  are  performed  in  program  AUTGEO,  whose  variables  are 
described  in  Table  1-3. 

2 

The  antenna  beam  patterns  all  follow  a  (sin  x/x)  shape  with  the 
beamwidth  as  a  parameter.  The  one-way  antenna  gain  relative  to  its  maxi¬ 
mum  value  is  given  by 

x  <  0.03 

x  >  0.03 


(  1  -  x2/6 

g(9)  -< 

I  (x  1  sin  x ) 2  , 
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Table  1-3.  Variables  in  Program  AUTGEO 


SYMBOL 

VARIABLE 

NAME 

DESCRIPTION 

ANTAZ 

Antenna  boresight  azimuth 

ANTEL 

Antenna  boresight  elevation 

bk 

BORE 

Unit  vector  along  antenna  boresight  axis 

pk 

RHO 

Vector  from  threat  to  target 

1  pk  1 

RANGE 

Range  to  target 

< 

AZTRUE 

True  target  azimuth 

K 

ELTRUE 

True  target  elevation 

RBMV 

Unit  vector  along  beam  axis,  S'E'Z'  system 

rk 

ABMV 

Unit  vector  along  beam  axis,  SEZ  system 

ak 

THETA 

Angle  of  target  off  beam  axis 

where  x  ■  1.4  9/9  and  the  beamwidth  is  29  .  When  9*9  ,  the  antenna 
o  oo 

gain  is  3  dB  down  from  its  value  at  9*0  or  g(9)  =  0.5.  Thus  to 
have  a  2°  beamwidth,  9q  *  1°  is  chosen. 

1.2.3  Target  (RPY)  Coordinate  System  and  Target  Geometry 

A  coordinate  system  attached  to  the  target  has  been  defined  to  specify 
the  locations  of  scatterers  and  emitters  attached  to  the  airframe  or  ship. 
An  aspect-dependent  cross-section  model  may  also  be  specified  in  this 
system. 

The  target  coordinate  system  has  its  origin  at  the  center  of  gravity 
(CG)  of  the  airframe  with  the  roll  axis,  R,  in  the  aft  direction,  the 

A  A 

pitch  axis,  P,  coward  the  right  side,  and  the  yaw  axis,  Y,  vertical 
when  the  target  is  level  relative  to  the  surface. 


The  location  of  any  point  in  the  target  system,  is  related  to 

its  location  in  the  SEZ  system  by  a  translation  and  a  rotation: 


pss  +  dee  +  Dzz 


“  PCG  +  PRR  +  PPP  +  PYY  * 


where  p„_  is  the  location  of  the  target  CG  in  the  SEZ  system. 

The  rotation  matrix,  M(0y,  6p,  0R)  of  the  coordinates  p^y  *n  the 
RPY  system  to  coordinates  PgEZ  in  the  SEZ  system  with  origin  at  the  tar¬ 
get  CG  is 
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The  translation  is  simply  pgEZ  »  pcG  +  pSEZ' 

The  angles  0R,  8p,  and  ©Y  are  the  roll,  pitch,  and  yaw  of  the 
target  in  the  SEZ  system.  The  velocity,  v,  of  the  target  is  aligned 
with  the  -R  axis  and  has  the  coordinates 

v  *  |v|(-cos  0^  cos  0p,  +sin  0Y  cos  0p>  sin  0p) 

The  velocity  vector  determines  the  orientation  of  the  target  coordinate 
system  except  for  the  roll  angle: 

sin  0p  -  v  •  Z  | v |  1 

tan  0y  =  v  •  E(-v  •  S)  1 

\ 

0R  »  °  . 

The  positive  sense  of  rotation  for  these  angles  is 

Roll  —  right  side  up, 

Pitch  —  forward  end  up, 

Yaw  —  right  turn. 
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Thus,  the  location  of  any  emitter  or  reflector  in  the  target 

at  time  k  is  given  in  the  SEZ  system  by 

rk  =  pCG  +  M(eR’eP’9Y)rk  * 

A  further  rotation  is  used  to  align  the  antenna  boresight  axis  with 
North  (-S) : 

(cos  E 

sin  Eq 

The  angles  Aq  and  Eq  are  the  azimuth  and  elevation  angles  of  the 
antenna  boresight  axis.  The  angles  of  the  source  at  rfc  projected  on 
the  azimuth  and  elevation  axes  are  then  easily  found  from 

tan  ejf  -  <r"  •  E) (-r£  •  S)'1 

sin  e£L  «  (r"-Z)  l^l'1  . 

The  range  difference  from  the  source  to  the  CG  along  the  line-of-sight 
from  the  antenna  is 

A  =  (rk  “  PCG^  '  PCG  ^PCgI 
1.3  RADAR  SIGNAL  PROCESSING 

The  digital  signal  processing  required  by  tracking  radars  employing 
pulsed  waveforms  is  described  in  this  section.  These  filters  are  used  in 
removing  noise,  detecting  the  amplitude  modulation  at  the  conical  scan 
frequency,  and  modeling  the  servo  motors  which  point  the  antenna.  The 
specific  descriptions  of  the  threat  models  and  ECM  controllers  of  Sections 
3  through  5  should  be  consulted  to  find  how  specific  filters  are  used 
within  a  model  and  what  their  parameter  values  are. 


* 


I 


f 


* 


t 


1.3.1  Digital  Filters 

Most  digital  filters  have  been  developed  from  transfer  functions 
written  for  analog  filters  in  the  continuous  frequency  domain  where  they 
are  described  by  their  gain  and  the  locations  of  their  poles  and  zeroes. 

The  transfer  function  for  the  analog  filter,  Ha(S),  is  transformed  into 
the  corresponding  discrete  filter  by  means  of  the  backward  difference 
operator,  where  T  is  the  sample  interval: 

TS  =  1  -  Z_1  . 

In  the  same  way  that  the  transfer  function  represents  a  differential  equa¬ 
tion,  the  resulting  system  function,  H(Z),  represents  a  difference  equa¬ 
tion  which  can  be  retrieved  to  provide  the  time-domain  representation  of 
the  filter  [3].  Table  1-4  lists  some  of  the  filters  used  in  the  simulation 
and  the  method  of  generating  the  output  sequence,  y^,  from  the  input 
sequence,  x,  . 


The  analog  transfer  function  has  the  form 


H  (S) 
a 


(2ojo/Q)(S+So) 

S2  +  (2u>o/Q)S+u£  ’ 


where  Q  =>  fQ(2Af)  *  and  fQ  »  ojo(2tt)  2  is  the  resonant  frqeuency  and 
Af  the  bandwidth  of  the  filter. 


The  same  transfer  function  represented  by  poles  and  zeroes  is 
Ha(S)  -  K(l  +  S/a)(l  +  S/b)“1(l+S/c)"1  , 

where 

b.c  -  ^[(f1  +<q"2-i)1/2] 

3  '  So 

K  -  2So(UoQ)-1  . 


» 
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Table  1-4.  Digital  Filters 


PROGRAM  ANALOG  DIFFERENCE  EQUATION 

POLEA  G(l+S/a)  ^  (l+aT)y,  =  y,  ,+G(aT)x, 

(ILPF1.LPF1)  K  k 

(ILPF2.LPF2) 

(ILPF3.LPF3 

POLEAA  G(l+S/G)-2  (l+aT)2yk  =  G(aT) 2xk+2(l+aT)yk_1 

"  yk-2 

FILZP  G(l+S/a) (1+S/b)-1  (aT)(l+bT)yk  =  G(bT) [ (1+aT) ^ 

-xk_1]+(aT)yk_1 

FILZPP  (l+S/b)(l+S/c)  (aT)(l+CT)(l+bT)yk=  (aT)(2+bT+cT)yk_1 

-(aT)yk_2+G(bT)(cT)[(l+aT)xk-xk_1] 

BPFZPP  2Ia/.^G(.S+-b.)  [l+(2aT/Q)+(aT)2]y,  =  2(l+aT/Q)y, 

s  +2Sa/Q+a  k  k~l 


+yk_2+2G(aT/Q) [(l+bT)xk-xk_1] 


Given  the  locations  of  the  poles,  b  and  c,  and  the  zero,  a,  the 
inverse  transformation  is 

2  , 

(i)  =  be 

o 

S  =  a 
o 

2u>  Q-1  =  b  +  c  =  Kbca-1  . 
o 

The  gain  is  not  an  independent  parameter,  so  that  the  transfer  function 
becomes,  upon  replacing  K; 

Hc(S)  =  (b+c) (S+a) (S+b)-1(S+c)_1  . 

The  discrete-time  version  of  the  filter  is  shown  in  Table  1-4. 

1.3.2  Automa tic  Gain  Control 


The  AGC  is  a  nonlinear  device  which  uses  a  linear  filter  to  produce 

its  output.  The  input  sequence,  z^,  is  averaged  by  a  linear  low-pass 

filter  to  produce  the  average  signal  z',  .  The  input  sequence  is  passed 

-1  * 

through  a  variable  gain  equal  to  z^  .  To  produce  the  AGC  output, 
z£  =  z^/z^,  as  shown  in  Figure  1-8. 

Programs  IAGCH  and  AGC4  implement  an  AGC  model  employing  a  low-pass 
filter  with  one  pole. 


Figure  1-8.  AGC  Block  Diagram 
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1.3.3  Antenna  Servomechanism  Models 


The  transfer  function  for  the  input  error  voltage,  v,  and  the  out¬ 
put  shaft  rotation  angle,  6,  is  assumed  to  be  of  the  form: 

S/v  =  K(l  +  S/u)ir1(l  +  S/ui3)  . 

The  corresponding  z-transform  obtained  from  the  replacement  TS  =  1-z  1 
is 


H(z)  =  K[l  +  (l-z“1)/(Tca1)]_1[l  +  (l-z“1)/(Tu,3)]  , 

which  corresponds  to  a  difference  equation  for  computing  the  output 
sequence,  5^,  given  by 

(Tu3)(l+Tu1)6k  =  KCTa)^  [(l+To)3)vk-vk_1]  +  (^3)6^  . 

The  shaft  rotations  are  integrated  by  a  filter  with  the  transfer 
function 


A(S) / 6 (S)  =  S  1  , 
whose  z-transform  is 

t, 

H(z)  =  T(l-  z-1)"1  . 

The  difference  equation  is  thus: 

\  3  Vi  +  T'\  • 

Programs  ISERV2  and  SERV02  implement  this  model. 

A  second  approach  to  implementing  the  servos  is  to  use  a  bilinear 
transform.  The  servo  motor  transfer  function  is  used  together  with  the 
bilinear  transform: 

TS  =  2(1-  z"1)(l  +  z"1)"1  . 
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to  obtain  the  corresponding  z -transform: 


-1 


H(z)  =  K 


(1-2/Toj3)  +  (1  +  2/Tu3)  6(z) 


z_1(l-  2/Ta^)  +  (1  +  2/To^)  v(Z> 


The  resulting  difference  equation  for  the  output  sequence  5^  is 

(l  +  2/Ta)1)6k+(l-2/Tw1)6k_1  =  K[(l  +  2/Tto3)vk+ (1- 2/Tw3)vk_1] 

Program  ISERV3  is  used  to  initialize  program  SERV02  when  this  servo 
model  is  desired. 


1.3. A  Error  Detector  in  Conical  Scan  System 

In  a  conical  scan  radar,  the  error  information  is  encoded  in  the 
amplitude  and  phase  of  the  modulation  of  the  detected  signal  at  the  conical 
scan  frequency.  This  -information  is  extracted  as  a  vector  error,  (o^.c^). 

The  phase  of  the  conical  system  is  defined  by 

(A^ef,E^ef)  =  {-sin(yk  +  <f>) ,  cos  (Yk  +  <P)  ]  , 

where  Yk  =  2tt  f ^  f^  k  is  the  true  phase  angle  of  the  conical  scanning  beam, 
<t>  is  a  phase  shift  to  account  for  the  phase  delay  in  the  error  signal  as 
it  is  processed  by  various  filters,  f  is  the  conical  scan  frequency,  and 
f  is  the  pulse  repetition  frequency. 

The  error  signals  for  each  component  of  the  angle  error  are  then 
obtained  from  the  error  signal,  Vk>  by: 


V  k+Z-1  , 

■>  ■  -t  & 


V  k+il-1  r 

e  -  -a  y  v"  Eref 

fck  £  L-l  C-l 

K  *  J-k  3  3 
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where  H  is  the  number  of  pulses  over  which  the  error  is  computed. 


Programs  IAEDET  and  AEDET  represent  this  model. 

1.4  THE  SCENARIO 

The  modules  related  to  the  scenario  to  be  exercised  are  described  in 
this  section.  The  target  models  are  given  in  Section  1.4.1.  These 
include  the  model  of  radar  cross-section,  the  ECM  antenna  model,  and  the 
ESM  antenna  model.  Section  1.4.2  describes  the  trajectories  of  the  tar¬ 
get  and  threat.  For  the  current  phase  of  the  project,  the  threat  is  sta¬ 
tionary.  Thus  only  the  target  trajectory  is  given.  The  ECM  objective  is 
described  in  Section  1.4.3.  This  objective  is  the  desired  angle  error  as 
a  function  of  time  which  the  ECM  system  is  attempting  to  induce  in  the 
threat  radars. 

The  models  describing  the  threat  portion  of  the  scenario  are  rele¬ 
gated  to  the  sections  describing  each  radar  type:  conical  scan  systems  in 
Section  2,  psuedo-monopulse  systems  in  Section  3,  and  monopulse  systems 
in  Section  4.  This  is  because  each  threat  model  is  specific  to  the  system, 
with  relatively  little  commonality. 

1.4.1  Target  Model 

The  target  model  includes  the  following  aspects  related  to  the  tar¬ 
get  platform: 

•  Radar  cross-section 

•  ESM  antenna 

•  ECM  antenna. 

Each  of  these  models  is  described  in  detail  below.  The  target  programs 
provide  an  interface  between  the  tracking  radar  and  the  programs  which 
model  the  ECM  techniques. 

In  this  phase  of  simulation  development,  the  weapon  employed  by  the 
threat  has  been  simplified  to  anti-aircraft  artillery  (AAA)  with  no 
pointing  angle  errors  of  its  own.  However,  the  simulation  is  capable  of 
incorporating  either  active  or  semi-active  homing  or  command-guided  mis¬ 
sile  weapons.  These  will  be  addressed  in  the  next  phase  of  the  project. 
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Radar  Cross-Section  Model 


The  target  radar  cross-section  can  be  chosen  from  one  of  three  avail¬ 
able  fluctuation  models:  deterministic,  exponential,  or  Rayleigh.  The 
correlation  time  of  the  fluctuations  can  also  be  controlled.  The  symbols 
in  the  program  RCSGEN  are  described  in  Table  1-5. 

The  distribution  of  cross-section,  c2>  when  the  exponential  model 
is  used,  is  given  by: 

_ i  -a2/°2 

f(02)  =  a2  •L  e  a2  >  0  . 

This  distribution  for  a2  is  generated  by  selecting  a  pseudo-random 
number,  r,  distributed  uniformly 

f  (r)  =  1.  ,  0  <  r  <  1  , 

and  then  making  the  transformation  C2  *  -c?2  £.nr.  It  has  the  average 
value  a2- 

The  distribution  of  cross-section,  a^,  when  the  Rayleigh  distribu¬ 
tion  is  selected,  is  given  by: 

— 2  -2a4^a4 

f(a^)  *  4a  4  a4  e  »  °4  0  • 

It  has  the  mean  a,  and  cumulative  distribution: 

4 

C  _  -2*/o4 

F(x)  =  J  C4  f(c4)dc4  «  1  -  (1  +  2x/c4)  e 
0 

The  distribution  for  c^  is  generated  by  selecting  a  pseudo-random  number, 
r,  distributed  uniformly  on  [0,1]  and  solving  the  equation  F(x)  “  r  for 
x.  This  is  done  iteratively  using  Newton's  method  in  the  program  RAYLGH. 


Table  1-5.  Variables  in  Radar  Cross-Section  Model, 
IRCSGN  and  RCSGEN 


SYMBOL 

VARIABLE 

NAME 

VALUE 

DESCRIPTION 

a 

o 

RCSDET 

0.0 

2 

Deterministic  cross-section,  m 

°i 

RCS2 

0.0 

2 

Average  exponential  cross-section,  m 

°4 

RCS4 

10.0 

2 

Average  Rayleigh  cross-section,  m 

T 

RCSTIM 

0.02 

Correlation  time  of  cross-section,  sec. 

f 

r 

PRFREQ 

400.0 

Pulse  repetition  frequency 

°2 

Cross-section  from  exponential  dis¬ 
tribution 

a4 

Cross-section  from  Rayleigh  distribu¬ 
tion 

°k 

Uncorrelated  cross-section  at  pulse  k 

°k 

Correlated  cross-section  at  pulse  k 

“I 

RCSMKV 

Correlation  time  filter  parameter 

♦ 


♦ 
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is  given  by  t  ■  kf^.  Thus,  the  parameter 
time  by 


The  monopulse  target  model  allows  multiple  point  sources  of  radar 
energy  which  can  be  either  reflectors  or  ECM  emitters.  The  jth  source 
is  located  in  the  airframe  RPY  coordinate  system  by  the  vector  r(j). 

The  cross-section  model  for  each  scatter  is  identical  to  that  for  point 
sources,  except  that  its  phase  may  be  specified  to  be  either  fixed  or 
chosen  uniformly  between  0  and  2tt.  For  those  sources  which  are  ECM 
emitters,  the  amplitude  and  phase  is  determined  by  the  ECM  control  algo¬ 
rithm. 

The  ECM  Antenna  Model 


where  the  time  at  sample  k 
a  is  related  to  correlation 


-<Tfr) 


-1 


The  ECM  antenna  model  interfaces  the  ECM  technique  controller  to  the 
tracking  radar  model  by  computing  the  emitted  power  in  appropriate  units 
and  allowing  some  flexibility  through  appropriate  input  parameters. 

The  directed  ECM  power  per  pulse,  emitted  by  the  ECM  antenna 

is  given  by 

ECM  .  ECM  .  gECM 
k  k  k  ’ 


where : 


pECM 

.ECM 


"k 


ECM  antenna  power  in  watts, 
ECM  antenna  gain, 

ECM  control  signal. 


The  current  model  for  G. 

k 

pointed  at  the  threat. 


is  a  constant,  modeling  an  antenna  always 


ESM  Antenna  Model 


The  ESM  model  includes  a  very  specific  model  for  the  power  per  pulse 
illuminating  the  target,  and  may  also  refer,  in  general,  to  any  other 
observables  which  are  used  by  specific  ECM  techniques.  Only  the  model 
for  a  passive  antenna  is  described  here. 

Given  that  power  density  p^  is  incident  on  the  ESM  antenna  at 
pulse  k,  the  observable  produced  by  the  antenna  is  the  power  received: 

_ESM  -ESM  aESM  „ 

Pk  =Gk  *  A  *  Pk  ’ 

ESM 

where  A  is  the  effective  antenna  aperture  and 
gain  of  the  antenna.  For  the  current  phase,  it  is 
tive  gain  of  the  ECM  antenna, 

_ESM  _ECM 
Gk  Gk  ’ 

modeling  ganged  antennas. 

Specific  ECM  techniques  should  be  consulted  for  a  description  of  any 
other  observables  used. 

1.4.2  Target  Trajectory 

The  current  interface  is  sufficiently  flexible  that  any  trajectory 
could  be  flown.  On  every  pulse,  an  entry  to  the  trajectory  generator 
is  made  to  compute  the  location  and  orientation  of  the  target.  The 
orientation  is  used  as  input  to  a  more  detailed  cross-section  model  only 
for  the  monopulse  threat. 

The  current  trajectory  generator  allows  a  straight-line  path  defined 
by  an  initial  position  and  velocity.  At  each  pulse,  k,  a  new  position 
is  computed  from 


FSM 

G,  is  the  relative 
k 

the  same  as  the  rela- 


t 


r 


k 


The  standard  trajectory  is  shown  in  Figure  1-9. 


♦ 


Figure  1-9.  Engagement  Scenario 
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1-32 


2. 


CONICAL  SCAN  SYSTEMS 


In  a  conical  scanning  radar,  there  are  one  or  more  beams  which  are 
either  mechanically  or  electrically  rotated  about  the  boresight  of  the 
antenna.  The  beam  rotation  causes  the  received  signal  to  be  modulated  at 
the  rotation  frequency  with  a  depth  of  modulation  which  grows  with  an 
increase  in  the  angle  between  the  target  and  the  antenna  boresight  and  a 
phase  related  to  the  direction  of  the  target  off  boresight.  The  resulting 
amplitude-modulated  signal  is  processed  to  retrieve  the  depth  of  modulation 
and  phase  as  input  to  the  servo  motors  which  move  the  antenna  boresight 
toward  the  target.  A  variant  of  the  basic  conical  scan  radar  is  conical 
scan  on  receive-only  (COSRO) ,  where  the  transmit  beam  is  fixed  at  the 

antenna  boresight  and  the  receiving  beam  performs  the  scanning.  This 

is  the  basic  system  considered  herein. 

This  section  describes  the  conical  scanning  threat  model  and  several 
ECM  techniques  which  have  been  developed  to  deny  tracking  to  the  threat. 

2.1  CONICAL  SCAN  THREAT  TRACKING  LOOP 

The  block  diagram  of  the  simulation  for  the  conical  scanning  threat 
model  is  shown  in  Figure  2-1.  The  major  signal  processing  functions  are 
shown  along  with  the  interface  between  the  threat  and  target  models. 

The  beam  geometry  requires  that  two  beams  be  defined,  one  for 
transmission  and  one  for  reception.  Figure  2-2  shows  the  beams  for  a 

CONSCAN  system  where  both  beams  coincide.  Figure  2-3  shows  the  beam 

geometry  for  a  COSRO  system.  The  latter  will  be  emphasized  in  this 
section. 

The  transmitting  beam  gain  in  the  direction  of  the  target  is  denoted 

T  T 

by  G  (ctk) .  The  power  density  illuminating  the  target  is  then 
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COSRO  Beam  Geometry 

where 


T 

P  =  transmitted  power  per  pulse. 


Pk  =  vector  from  threat  to  target. 


The  signal  due  to  the  skin  return  and  the  ECM  energy  emitted  which 
is  detected  by  the  receiving  beam  is 


2k  = 


,  „  n !/2  ^  dECM 

(pk°k>  +  Pk 


1/2 


i*t> 


AR  GR(aR)[GR(0)]_1(4TT 


lpk!2)'1+wk 


(2.2) 


where 


ak  =  target  cross-section, 

ECM 

=  directed  ECM  power  per  pulse 

£ 

A  =  effective  perture  of  receiving  beam, 
£ 

G  =  gain  of  receiving  beam. 


The  angle  4>  is  randomly  chosen  from  a  uniform  distribution  between  0 
and  2tt.  The  additive  noise,  wk,  is  Rayleigh-distributed: 


f(w  )  =  4wk  w  2  exp{-2wk/w}  .  (2.3) 

i 

The  detected  signal,  zk>  then  enters  the  AGC,  the  error  signal 
filter,  the  error  signal  detector,  and  the  servo  motors.  All  of  these 
signal  processing  models  are  fully  described  in  Section  1.3.  The  parameter 
values  for  the  standard  conical  scan  threat  are  shown  in  Table  2-1. 


2.2  SIMULATION  OF  ECM  TECHNIQUES 

The  expressions  for  the  required  ECM  emitter  control  have  been 
developed  in  a  previous  report  [2].  In  order  to  implement  them,  means 
for  computing  the  required  signals  must  be  developed,  in  particular  for 
the  average  values  of  modulated  signals. 


Table  2-1.  Standard  COSRO  Threat  Parameters 


THREAT  COMPONENT 

AND  PROGRAM  NAME 

PARAMETER 

VALUE 

SYMBOL 

NAME 

Conical  Scan 

Frequency 

25  Hz 

f 

c 

CSFREQ 

ICSG,  CCS 

Phase  Shift 

-12°. 70 

Transmitter 

p.r. f . 

400  Hz 

fr 

PRFREQ 

THRXMT 

Power 

200  kW 

PT 

PTRANS 

(Beam  #2) 

c  .  „  (  COSRO 
Squint  ; 

( CONS CAN 

o°.o 

0° .  6 

8 

BETA 

Phase 

0°. 

yD 

GAMMA 

Gain 

40  dB 

gt(0) 

RELGN 

Half  | COSRO 

Beamwid  th  | CON S  CAN 

1°  .5 

Is.  0 

0o 

BMWID 

Reciever 

Squint 

0° .  6 

6 

BETA 

THRRCV 

Phase 

o 

o 

O 

vo 

GAMMA 

( Beam  01 ) 

Half 

Beamwid th 

1°.0 

So 

BMW  ID 

Aperture 

0.88  m2 

ar 

RCVARE 

Noise 

-130.0  dB 

W 

AGC 

Gain 

1.0 

AGC  5 

Pole 

5.0  Hz 

Error  Signal  Filter 

Zero 

10  Hz 

ERFIL 

Resonant  Freq. 

23  Hz 

Q 

0.92 

Error  Signal  Detector 

Averaging  Period 

1  Con.  Scan 

Prd. 

AEDET 

Gain 

0.046 

Servos 

ISERV02 ,  SERV02 

Zero 

1.6  Hz 

Pole 

0.16  Hz 

Gain 


2-6 


200.0 


The  following  sections  describe  the  computer  programs  which  are  used 
to  simulate  ECM  controllers  against  COSRO.  The  ECM  technique  is  divided 
into  two  functions:  (1)  the  generation  of  a  waveform,  and  (2)  control  of 
the  waveform  by  a  control  algorithm. 

2.2.1  Parametric  Controller  Against  COSRO 

The  minimum-power  ECM  signal  required  to  induce  a  constant  received 
power  at  the  tracking  radar  has  been  developed  in  [2]  and  is  given  by: 

(2. A) 

The  simulation  of  this  ECM  technique  consists  of  three  parts:  part 
one  is  the  computation  of  the  control  voltage  using  the  estimates  for  the 
unknown  quantities,  part  two  develops  estimates  for  the  required  signals, 
and  part  three  is  an  adaptive  scheme  for  estimating  unknown  threat  para¬ 
meters.  Each  part  is  described  below. 

Because  the  COSRO  transmitting  beam  does  not  rotate,  the  ESM  power 
received  per  pulse  does  not  contain  as  much  information  as  the  modulated 
CONSCAN  observable.  It  will  be  assumed  that  there  is  available  an  addi¬ 
tional  observable,  the  true  total  tracking  angle  error  at  the  time  of 
each  pulse  transmission.  The  observables  are  thus: 

Y^  *  powt-’L'  received  by  ESM  antenna  on  pulse  k, 
y^  *  tracking  angle  error  on  pulse  k,  . 

The  variables  used  in  the  computer  program  are  described  in  Table  2-2. 
Some  of  the  parameters  are  important  properties  of  the  threat  radar,  such 

A  A  A 

as  Jo,  S,  and  0  ,  which  are  not  determined  adaptively  by  the  estima¬ 
tion  scheme.  Sensitivity  of  the  techniques  to  the  assumed  values  is  dis¬ 
cussed  in  the  performance  evaluation. 


\ 


\  Eltk1 

E<V 


•1+nax  [E{GR(o£))]/E{GR(ciR)} 
k 
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Table  2-2.  Variables  in  Programs  CSR01E.F4  and  CSR02E.F4 


SYMBOL  VARIABLE  DESCRIPTION  VALUE 


/s 

?. 

LPULSE 

Number  of  pulses  per  conical 
scan  period 

16 

✓V 

6 

EBETA 

Estimated  squint  angle 

0° .  6 

/\ 

9 

0 

EBWID 

Estimated  receiving  half 
beamwidth 

H* 

o 

O 

Nd 

NSETTL 

Number  of  pulses  to  delay 
start  of  ECM 

100 

a 

RISE 

Angle  error  filter  constant 

0.75 

0 

FMAG 

Final  desired  angle  error 

0.015 

rad 

> 

a 

AFILT 

ESM  power  filter  constant 

0.8 

g 

GAIN 

Estimation  technique  parameter 

0.0 

Initialization 

SYMBOL 

VARIABLE 

DESCRIPTION 

VALUE 

\ 

APLAT 

Target  platform  constant 

A-1  " 

500 

9k 

ADES 

Desired  angle  error  at  time  k 

9-l  " 

0.0 

k 

KPULSE 

Pulse  counter 

-1 

Pk 

KP 

Conical  scan  counter 

1 

Yk 

CSANG 

Conical  scan  phase  angle 

0.0 

d 

Y 

DCSANG 

Increment  to  per  pulse 

2UJT1 

Sk 

SD 

Standard  deviation  of  estimate 
for  Ak 

S-2  * 

S-1 

\ 

YANG 

Average  received  ESM  power 

Y-1  * 

0.0 

Yk 

ANG 

True  tracking  angle  error 

9-l  * 

0.0 
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COSRO  Controller 


The  controller  computes  the  power,  u^,  to  be  applied  to  the  ECM 
emitter  on  the  kth  pulse  using  Eq.  2.4  as  a  basis.  The  estimates  used 
for  the  various  quantities  are: 


E[Y^3/E[Yk]  -  Yk 

max  {E[GR(aR)]}  ■+•  G(a^in) 

E[GR(a£)]  -  G(ak)  . 

A 

An  estimate  for  the  target  platform  parameter,  Ak>  has  been  used 
instead  of  its  true  value  at  time  k.A^.  The  estimate  may  be  either  a 
constant  chosen  before  the  engagement  begins  and  which  does  not  change 
during  the  engagement  (non-adaptive)  or  it  maybe  determined  during  the 
engagement  by  an  adaptive  estimator.  Both  techniques  have  been  imple¬ 
mented.  Evaluations  are  given  in  Section  5. 

The  use  of  Yk>  the  average  ESM  power  received  per  pulse,  is  valid 
in  this  case  because  Yk  is  slowly  varying  in  time,  since  the  COSRO 
transmitting  beam  is  not  squinted  off  boresight. 


The  maximum  receiver  gain  occurs  when  the  angle  between  the  target 

and  beam  axis,  a.  ,  is  a  minimum.  The  superscript  'R'  is  dropped  because 

inin 

all  references  will  be  to  the  receiving  beam.  The  angle  is  defined 

by 

ml  n  R 

c^ln=min  .  (2.5) 


Thus,  the  ECM  power  transmitted  is  given  by 


“k = 


Yk[-l  +  G(ot“ln)/G(ak)]  , 


k  <  N , 
a 

otherwise 


(2.6) 


in  the  computer  model. 

The  technique  employs  a  feedback  controller  because  the  observations, 
Yk>  of  the  radar  power  illuminating  the  target  are  used  to  adjust  the 
strength  of  the  ECM  signal,  i^.  The  technique  becomes  adaptive  when  an 

A 

estimation  scheme  is  used  to  determine  the  value  of  ^  while  the  tracked 
radar  is  being  controlled. 

The  filter  for  Yk  runs  for  pulses  in  order  to  obtain  a  stable 

estimate  during  which  no  ECM  power  is  transmitted. 

COSRO  Signal  Generators 

Each  time  the  target  is  illuminated  by  a  pulse  from  the  tracking 
radar,  the  internal  conical  scan  phase  angle  is  updated  by: 

Yk  *  Yk-1  +  dY  *  (2,7) 

and  the  average  value  for  ESM  power  received  is  updated: 

\  “  a'  Yk-1  +  U-a>  \  •  (2.8) 

That  Yk  does  not  correspond  to  the  "true"  phase  angle  in  the  tracking 
radar  is  not  important.  It  is  used  only  to  generate  an  ECM  signal  of  the 
current  frequency  with  the  proper  amplitude  variation. 

Once  per  control  scan  period,  a  new  value  for  the  desired  angle 
error,  @k,  is  obtained  from  the  previous  value  and  the  final  desired 
angle  error,  9: 
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(2.9) 


0^  »  ot  9k_J  +  (1  -  a; j  . 


The  effect  is  to  gradually  increase  the  desired  angle  error  from  zero  to 
its  final  value. 

The  internal  version  of  the  angle  between  the  target  and  the  receiving 
beam  axis  is  given,  for  small  angles,  by: 


\  -  i  *  »2  -  s  “■  \ 


|ek  -  3 |  . 


(2.10) 


The  beam  shape  G  is  then  used  to  compute  G(t\)  anc*  G^^111),  where 


G(a)  *  (x  1  sin  x)2  , 


x  =  1.  A  a  0  1  , 
o 


(2.11) 


and  S  is  the  estimated  half  beamwidth. 
o 


Once  the  estimate,  A^,  is  obtained,  it  is  a  simple  matter  to 
compute  the  ECM  signal,  u^. 

Adaptive  Estimate  for  Platform  Parameter 

A  new  parameter  estimate,  A^,  is  determined  only  once  per  conical 
scan  period,  when  the  value  of  k^  »  1+k  mod  £  =  1.  That  estimate  is 
then  used  throughout  the  next  conical  scan  period  to  generate  the  ECM 
signal.  When  k^  =  1  and  k  21 


Pk  =  Pk-1  +  1  * 


(2.12) 


and  a  new  estimate,  A^,  is  computed  by  comparing  the  desired  angle 
error,  9^,  with  the  actual  angle  error,  y^: 
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1 


\  =  ‘Sc-£  "  P  Sk-2£ 
k-1 


l(x'  (l"v)’ 


(2.13) 


»here!  sk  • 


Parametric  Controller:  Non- Adaptive  ECM  Technique  Against  COSRO 

This  ECM  technique  is  very  similar  to  the  adaptive  technique  described 
above  but  no  attempt  is  made  to  ascertain  the  value  of  the  platform  para¬ 


meter,  A^,  during  the  encounter.  Instead,  a  pre-specif ied  value  of 
A^=Ao  is  used  throughout  the  engagement.  The  controller  still  employs 
feedback  through  the  use  of  its  observations,  and  y^  to  adjust  for 

illumination  by  the  tracking  radar. 

The  input  parameters  are  the  same  as  in  Table  2-2  except  that 
A^=Aq  =  500.0  is  also  included.  Since  no  estimation  technique  is 
employed,  is  not  used.  The  initialization  is  otherwise  the  same 

as  in  Table  2-2. 

2.2.2  Input/Output  Controller  Interface  for  COSRO 

Two  techniques  were  used  to  interface  the  jamming  module  with  controllers 
designed  for  linear  systems.  Both  techniques  use  the  actual  radar  boresight 
angular  error  (with  additive  gaussian  noise  if  desired)  as  the  state  parameter 
y.  As  each  new  value  of  y  is  presented  to  a  controller  (or  identifier- 
controller  pair)  the  controller  computes  a  new  control  u  which  is  used  to 
control  the  ECM  signal  transmitted  during  the  next  sampling  interval. 

A  reference  signal  is  used  which  is  modulated  at  the  radar's  conical 
scan  frequency.  Both  of  the  interface  techniques  described  here  use  sinusoidal 
modulation  at  exactly  the  scan  frequency.  The  assumption  is  that  the  scan 
rate  is  fixed  and  known.  For  the  "depth  modulation"  technique  the  power  to 
be  transmitted  as  ECM  during  pulse  k  is  given  as  a  function  of  u: 


ECMPOW,  -  1  -  u  sin (2uk/ £)  , 
k 


where  2,  is  the  number  of  pulses  per  conical  scan.  For  the  "amplitude 
modulation"  technique: 

ECMPOWk  =  |u|  -  u  sin(2mk/2)  . 


Sampling 

ECM  power  is  generally  transmitted  on  each  pulse  thus  requiring  one 
of  the  above  calculations  to  be  performed  for  each  pulse.  However,  it  is 
undesirable  to  sample  the  angle  error  y  for  use  in  identification  and 
computation  of  u  (the  control  parameter)  at  a  rate  faster  than  the 
dynamics  of  the  radar  system.  Thus  some  subsampling  is  required.  It  is 
also  apparent  from  the  modulation  equations  that  u  should  be  held  constant 

A 

for  a  multiple  of  2.  pulses  (integer  number  of  conical  scans).  In  general 
with  the  ARMA  based  controllers  which  have  been  implemented  against  the 

A 

"standard"  threat  a  sampling  interval  of  22,  pulses  has  been  satisfactory. 
Ramping 

It  has  been  found  that  when  control  is  first  applied  to  the  radar  it  is 
advantageous  to  preempt  drastic  control  commands  which  can  cause  the  radar  to 
lose  track  or  behave  erratically.  This  is  particularly  true  if  no  identifi¬ 
cation  phase  is  used  or  if  this  phase  is  short.  For  identification/control 
mechanisms  this  is  easily  accomplished  by  smoothly  adjusting  Y  (desired 

boresight  angle  error)  from  zero  (identification  phase)  to  its  final  value. 

An  exponential  ramp  with  a  time  constant  of  approximately  1  second  is  generally 
adequate. 

For  certain  direct  controllers  such  as  output  matching,  no  identification 
phase  is  involved  which  means  that  no  estimate  of  the  closed  loop  gain  of  the 
system  is  available  to  startup.  In  this  case,  it  is  desirable  to  limit  the 
control  in  the  beginning.  Note  also  that  since  the  depth  modulation  tech¬ 
nique  transmits  a  constant  average  power  (typically  10  dB  or  so  above  the  skin 
return)  it  is  desirable  to  smoothly  turn  the  ECM  power  up  when  using  this 
technique.  These  goals  are  easily  accomplished  by  limiting  the  ECM  power 
transmitted  to  the  level  of  an  exponential  ramp  as  used  above  for  YREp. 
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3.  PSUEDO-MONOPULSE  SYSTEMS 


A  psuedo-monopulse  system  consists  of  two  radar  beams,  n  radians 
out  of  phase,  each  squinted  off  the  antenna  axis  by  angle  2,  which  are 
either  electrically  or  mechanically  rotated  about  the  antenna  axis.  One 
of  the  beams  transmits  and  both  of  them  receive.  The  beam  configuration 
is  shown  in  Figure  3-1.  The  following  sections  describe  the  standard  threat 
model  and  the  ECM  techniques  which  have  been  developed. 

3.1  PSUEDO-MONOPULSE  TRACKING  LOOP 

The  block  diagram  of  the  simulation  for  psuedo-monopulse  threat  model 
is  shown  in  Figure  3-2.  The  major  signal  processing  functions  are  shown 
along  with  the  interface  between  the  threat  and  target  models. 

The  beam  geometry,  shown  in  Figure  3-1,  requires  that  two  beams  be 
defined.  One  is  the  main  beam,  which  both  transmits  and  receives.  The 
other  is  the  auxiliary  beam,  which  only  receives. 

The  transmitted  power  density  illuminating  the  target  is  given  by 

pk  =  PT  GT(cx£)  (4tt  |pk|2)'1  ,  (3.1) 

T 

where:  P  =  transmitted  power  per  pulse, 

T 

G  =  main  beam  gain, 

T 

ak  ■  angle  of  target  from  main  beam  axis, 

Pk  =  range  to  target. 

The  signals  detected  by  each  beam  have  contributions  from  both  the 
power  reflected  from  the  target  and  from  ECM  power: 
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elevation  axis 

MAIN  BEAM 
3  dB  CONTOUR 


Figure  3-1.  Psuedo-monopulse  Beam  Geometry 


3-2 


monopulse  Threat  Model  Block  Diagram  Showing  Program  Names 


where : 


1/2  2 

(pkak)1/2  +  PkCM  aBgB(\}  [GB(0)  ]-1(4tt  |pk|V1  +  wk  , 


ak  =  target  cross-section. 


Pk  =  directed  ECM  power  per  pulse, 

B 

A  =  effective  aperture, 

g 

G  =  gain  of  beam, 

g 

ak  =  angle  of  target  from  beam  axis, 
B  =  main  (T)  or  auxiliary  (R)  beam. 


The  angle  <j>  is  randomly  chosen  from  a  uniform  distribution  on  the  inter¬ 
val  [ 0 , 2tt ] .  The  additive  noise,  wk>  is  chosen  from  a  Rayleigh  distri¬ 
bution  (see  Eq.  2.3). 

Each  of  the  detected  signals  enters  an  AGC.  The  difference  of  the  AGC 
outputs  enters  an  error  signal  filter,  the  error  signal  detector,  and  the 
servo  motors.  All  of  these  signal  processing  models  are  fully  described 
in  Section  1.3.  The  parameter  values  for  the  standard  psuedo-monopulse  threat 
model  are  shown  in  Table  3-1. 

3.2  SIMULATION  OF  ECM  TECHNIQUES 

The  expressions  for  the  required  ECM  emitter  control  have  been 
developed  and  analyzed  in  [2j.  The  following  sections  describe  the  com¬ 
puter  programs  which  are  used  to  simulate  ECM  controllers  against  psuedo-monopulse. 
Again  each  ECM  technique  is  divided  into:  (1)  the  generator  of  a  waveform, 
and  (2)  the  control  of  the  waveform  parameters. 

3.2.1  Parametric  Controller  Against  Psuedo-monopulse 


The  control  is  the  voltage,  uk>  of  the  ECM  emitter  on  the  kth  pulse. 
The  periodic  waveform  is: 
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Table  3-1.  Standard  Psuedo -monopulse  Threat  Model  Parameters 


THREAT  COMPONENT 

AND  PROGRAM  NAME 

PARAMETER 

VALUE 

SYMBOL 

Conical  Scan 

Frequency 

25  Hz 

f 

c 

ICSG,  CSG 

Phase  Shift 

-13°. 29 

^REF 

Transmitter 

p. r. f . 

400  Hz 

fr 

THRXMT 

Power/Pulse 

200  kW 

PT 

(Beam  #2) 

Squint  angle 

0° .  6 

S' 

Phase 

180° 

Y0 

Gain 

40  dB 

GT(0) 

Half  Beamwidth 

1°.0 

9o 

Aperture 

0.88  m2 

Receiver 

Squint  angle 

o 

0 

o> 

6 

THRRCV 

Phase 

O 

o 

O 

Vo 

Half  beamwidth 

1°.0 

0o 

Aperture 

0.88  m2 

Noise 

-130.0  dB 

w 

AGC 

Gain 

1.0 

AGC4,  AGC5 

Pole 

0.5  Hz 

Error  Signal  Filter 

Zero 

10.  Hz 

ERFIL 

Resonant  freq. 

23.  Hz 

q 

0.92 

Error  Signal  Detector 

Averaging  Period 

1  Conical 

Scan  Period 

AEDET 

Gain 

0.023 

Servos 

Zero 

1.6  Hz 

1SERV2,  SERV02 

Pole 

0.16  Hz 

Gain 

200.0 

NAME 


CSFREQ 

PRFREQ 

PTRANS 

BETA 

GAMMA 

RELGN 

BMWID 


BETA 

GAMMA 

BMWID 

ANOISE 
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w 


s  0  ’  k  1  Nd 

| 6(k  mod  i)  ,  otherwise 


(3.3) 


where  is  a  delay  before  the  ECM  is  turned  on  and  is  the  emitted 

amplitude.  The  waveform  has  period  i. 


The  amplitude  of  the  emitted  pulse,  S^,  is  controlled  by  the  out¬ 
put  of  a  filter  whose  input,  y  ,  is  the  true  observed  tracking  angle 
error  on  the  kth  pulse: 


yk  =  “s  yk-l  +  (1_as)yk 


The  amplitude  is  feedback-controlled  according  to 


— 2fl~2 

=  SECM  e-yk  o 
k 


(3.4) 


(3.5) 


ECM 

where  S  is  a  constant  and  6q  is  a  parameter  representing  the  half- 
beamwidth  of  the  tracking  radar.  The  parameters  for  the  parametric 
controller  against  psuedo-monopulse  are  shown  in  Table  3~2. 


3.2.2  Input/Output  Controller  Interfaces 


As  in  the  COSRO  simulation  the  observable  radar  state  parameter  y  is 
taken  to  be  the  actual  boresight  angular  error  of  the  pseudo-monopulse  radar. 
A  simple  one  pole  smoothing  filter  is  used  with  a  time  constant  of  200  pulses 
to  smooth  the  rather  irregular  behavior  of  this  parameter.  This  smoothed 
estimate  is  sampled  and  fed  to  the  desired  identification/control  setup  which 
then  computes  a  new  control  u  to  be  used  during  the  next  sampling  interval. 
The  sampling  interval  chosen  for  this  process  (see  Section  2.2.2)  is  120 
pulses/ sample . 


The  actual  waveform  transmitted  as  ECM  is  defined  in  Section  3.2.1, 
Eq.  (3.3).  The  pulse  amplitude  S^  of  the  pulse  is  determined  from  the 
control  value  u  by: 
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Table  3-2.  Description  of  Variables  in  Parameteric  Controller 
Against  Psuedo-Monopulse 


SYMBOL 

VARIABLE 

VALUE 

DESCRIPTION 

a 

ALFCOM 

0.0 

Comb  filter  constant 

Nd 

NSETTL 

100 

Number  of  pulses  to  delay  start  of  ECM 

i 

LHAT 

16 

Estimated  number  of  pulses  in 
conical  scan  period 

9 

o 

EBWID 

Estimated  half  beamwidth 

sECM 

SECM 

1000.0 

ECM  voltage 

i 

IECM 

9 

ECM  pulse  spacing 

a 

s 

SFILT 

0.995 

Angle  error  filter  parameter 

k 

KPULSE 

-1 

Pulse  counter 

\ 

SANGER 

0.0 

Filtered  angle  error 

KMODL 

Estimated  conical  scan  pulse  position 

ki 

KMODE 

Pulse  number  relative  to  ECM  pulse 
spacing 

\ 

YK(2) 

Observed  true  tracking  angle  error 

yk 


YK(1) 


Observed  ESM  power  received 


(U  +  UBIAS)(1---"  >  * 


The  second  term  in  the  above  product  is  simply  an  exponential  ramp  as 
used  against  COSRO.  The  term  is  an  estimate  of  the  "critical  u" 

at  which  the  radar  is  just  becoming  unstable.  This  term  is  not  really 
necessary  for  certain  controllers  which  can  determine  the  bias  term  them¬ 
selves  . 


4. 


MONOPULSE  SYSTEMS 


A  monopulse  system  has  at  least  two  feed  horns  in  the  tracking  direc¬ 
tion  which  are  not  rotated,  but  rather  sense  the  wavefront  distortion 
coherently,  resulting  in  an  error  signal.  The  particular  monopulse  system 
simulated  is  a  full  three-channel  amplitude  comparison  monopulse. 

This  section  describes  the  monopulse  threat  model  and  the  ECM  tech¬ 
niques  which  have  been  developed  to  counter  this  threat. 

4.1  MONOPULSE  THREAT  TRACKING  LOOP 

The  monopulse  threat  model  is  intended  to  represent  the  behavior 
or  a  true,  three-channel  monopulse  antenna,  although  the  parameters  used 
do  not  represent  a  particular  system.  The  monopulse  model  differs  from 
the  conical  scanning  systems,  with  or  without  compensation,  in  that  a 
spatially  extended  target  is  required.  The  associated  geometry  has  been 
described  in  Section  1.2.3. 


The  block  diagram  for  simulation  of  the  monopul se  radar  is  shown 
in  Figure  4-1,  which  includes  the  program  names  and  the  interface  to  the 
target  model. 

There  are  four  beams  defined  for  the  antenna,  each  squinted  off-axis 
by  the  same  amount  and  located  at  phase  angles  corresponding  to  the  azimuth 
and  elevation  axes,  as  shown  in  Figure  4-2. 


The  power  density  at  the  target  transmitted  by  the  tracking  radar  is 


Pk  =  P 


(Att  IpJ2)"1 


(4.1) 
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Figure  4-1.  Monopulse  Threat  Model  Block  Diagram  Showing  Program  Names 


Figure  4-2  Monopulse  Beam  Geometry 
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where : 


power  transmitted  per  pulse. 


=  gain  of  jth  beam, 

o^j)  =  angle  of  target  from  jth  beam  axis, 
Pk  =  vector  from  threat  to  target. 


An  extended  target  is  modeled  as  a  set  of  point  sources  (scatterers 
or  emitters),  each  characterized  by  an  amplitude  and  phase 

The  amplitude  is 


1/2 

[ak(i)  •  pk  •  g(ak(i))]  scatterer 

1/2 

[uk(i)  •  g(CL  (i))]  emitter 


(4.2) 


where  g  is 
on-axis ,  and 
The  phase  of 
scatterers , 

[  0 , 2tt  ] , 


the  relative  gain  of  the  receiving  beam  normalized  to  1 
<*k(i)  is  the  angle  of  the  source  off  antenna  boresight. 
the  ECM  emitters  is  determined  by  the  controller.  For 

can  be  either  fixed  or  chosen  uniformly  on  the  interval 


The  signal  received  in  the  Rth  channel  is 


where 


S£  =  A(4tt  |pk|2)_1  E  fja^  cos[4>i-(|>j  -  27rK(Ai  -  A^)] +wfc  ,  (4.3) 

1  Z  =  sum  channel 

A  ^ 

yk  (i)  Z  *  azimuth  channel 

EL 

yk  (i)  Jc  »  elevation  channel, 

A  *  antenna  effective  aperture, 

A^  *  range  to  ith  source, 

K  *  A  \  X  is  transmitted  wavelength. 


The  noise,  wk>  is  Rayleigh  distributed  (see  Eq.  2.3). 
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The  sum-channel  signal  enters  a  low-pass  filter  whose  output  is  used 
to  control  the  gain  of  the  azimuth  and  elevation  channels.  Those  gain- 
controlled  signals  are  then  filtered  and  used  to  drive  to  the  servo  motors. 

The  parameters  of  the  standard  monopulse  threat  model  are  shown  in 
Table  4-1. 

4.2  SIMULATION  OF  ECM  TECHNIQUES 
Input/Output  Interface  for  Monopulse 

The  observable  parameter  for  the  monopulse  radar  simulation  is  the 
azimuth  component  of  the  pointing  error.  This  is  in  contrast  to  the  COSRO 
and  pseudo -monopulse  simulations  in  which  the  observable  (y)  is  the  total 
angular  error.  This  value  is  subsampled  by  a  factor  of  30  and  used  to  compute 
a  control  parameter  u  by  the  identification/ control  scheme  selected.  This 
value  is  then  used  to  compute  the  amplitude  and  phase  of  ECM  transmitted  from 
each  of  2  antennas  during  the  next  30  pulses. 

First  the  control  is  "linearised"  by  defining  a  gain  term  A: 


i  i 

u  if  |  ]  u  1  —1 j  <  . 4 
1.4SGN(u)  otherwise. 

A  phase  is  also  computed: 

|  it  if  |u  |  >  1 . 

(j)  *  2-rr  •  WAVNUM  •  AR  + 

( 0  otherwise, 

where  WAVNUM  is  the  wave  number  of  the  output  and  AR  is  the  difference 
in  range  of  the  radar  from  the  two  ECM  antennas.  This  equation  means  that 
the  output  from  the  two  antennas  will  arrive  at  the  radar  with  a  phase 
difference  of  t  if  |u|  >  1.  and  if  iu!  <.  1. 

Now  the  amplitude  and  phase  of  the  output  can  be  written: 

AMP(l)  *  1. 

AMP (2)  -  A  •  AMP(l) 

PHASE(l)  -  0. 

PHASE (2)  -  $ 
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u*  +  l 


lu*  -  1 


where  u 


*  J 


Table  4-1.  Standard  Monopulse  Threat  Model  Parameters 


THREAT  COMPONENT 

AND  PROGRAM  NAME 

PARAMETER 

VALUF 

SYMBOL 

NAME 

Transmitter 

p.r.f . 

400  Hz 

f 

r 

PRFREQ 

THRXKT 

Power/ pulse 

200  kW 

pT 

PTRANS 

Squint  angle 

0° .  6 

3 

BETA 

Phase 

0°, 90° , 180° , 
270° 

Y. 

u 

GAMMA 

Gain 

40  dB 

GNAXIS 

Half  beamwidth 

1°.0 

e 

0 

BMWID 

Reciever 

Noise 

-130  dB 

w 

ANOISE 

THRRCV 

A-l 

10  m"1 

K 

WVNMBR 

Aperture 

0.88  m2 

A 

RCVARE 

Sum  Channel 

Gain 

10.0 

ILPF1,  LPFl 

Freq.  Cutoff 

2.0  Hz 

Azimuth  Channel 

Gain 

10.0 

1LPF2,  LPF2 

Freq.  Cutoff 

2.0  Hz 

Elevation  Channel 

Gain 

10.0 

ILPF3,  LPF3 

Freq.  Cutoff 

2.0  Hz 

Servos 

Zero 

1.6  Hz 

Pole 

0.16  Hz 

Gain 

200.0 

5.  SAMPLE  RUNS 


The  simulation  of  the  adaptive  closed-loop  ECM  techniques  described  herein 
accomplished  the  primary  objective  of  algorithm  development.  That  is,  the 
algorithms  described  in  [2]  and  this  report  were  developed  in  part  by  using 
the  simulation  to  understand  ECM  behavior.  There  was  never  any  intention  to 
use  the  current  simulation  for  "production  runs".  Early  sample  runs  of  the 
simulation  are  meaningless  for  technique  evaluation.  Therefore  the  limited 
number  of  sample  runs  given  in  this  section  utilize  the  algorithms  in  their 
final  forms  (as  described  in  Sections  2  through  4). 

As  the  controllers  will  be  developed  further  in  subsequent  phases  of  the 
program,  even  those  runs  presented  herein  are  not  indicative  of  the  ultimate 
gains  to  be  achieved  using  adaptive  closed-loop  control  techniques.  For  this 
reason  only  a  sampling  of  the  runs  is  presented.  They  are  meant  to  be 
representative  of  the  performance  of  the  techniques  at  the  current  time,  and 
contain  a  variety  of  parameters,  methods,  and  tests.  Definitive  conclusions 
*■  from  such  runs  generated  by  the  simulation  will  be  given  during  the  next  phase 

of  the  program,  when  algorithms  will  be  developed  to  work  against  actual  threat 
radar  systems. 

I  Sample  runs  are  given  against  all  three  radar  systems:  COSRO  in  Section 

5.1,  psuedo-monopulse  in  Section  5.2,  and  monopulse  in  Section  5.3.  Parametric 
controller  results  are  shown  against  COSRO  and  psuedo-monopulse,  and  input/output 
model  controllers  are  shown  against  COSRO  and  monopulse.  These  choices  are 
♦  due  to  the  following  considerations.  COSRO  is  the  simplest  system,  due  to 

its  direct  coupling  of  observation  and  ECM  and  its  nearly  linear  structure. 

Thus  parametric  and  input/output  controllers  were  developed  against  COSRO. 

The  emphasis  on  COSRO  is  not  so  much  indicative  of  threat  system  importance, 

,  but  rather  shows  the  "walk  before  you  run"  research  performed  under  the  current 

project.  Then  a  significant  amount  of  analysis  (see  [2])  was  performed  to 
understand  the  technique  used  against  psuedo-monopulse,  and  only  later  was  the 
control  tried.  The  first  controller  tried  with  the  new  technique  was  the 
parametric  controller.  Meanwhile,  the  input/output  controller  work  emphasized 
monopulse. 
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In  all  cases  only  brief  descriptions  are  given  of  the  controllers,  and 
there  is  no  attempt  at  global  conclusions,  due  to  the  developmental  nature 
of  the  algorithms.  In  all  cases  the  controllers  "worked",  in  the  sense  that 
they  induced  the  error  behavior  desired.  This  control  of  the  threat  radar 
indicates  the  power  in  the  approach. 

5.1  SAMPLE  RUNS  AGAINST  COSRO 

Since  most  of  the  algorithm  development  was  done  for  the  COSRO  radar, 
the  most  sample  runs  are  presented.  Table  5-1  describes  the  sample  runs, 
with  Figures  5-1  through  5-20  giving  outputs.  In  this  and  subsequent  tables 
the  plot  types  are  given  by: 

1  -  Angle  error  versus  time  at  start  of  ECM 

2  -  Angle  error  histogram  at  start  of  ECM 

3  -  Angle  error  versus  time  in  steady-state 

4  -  Angle  error  versus  histogram  in  steady-state 

5  -  ECM  signal  versus  time  (envelope) 

6  -  ECM  histogram 

7  -  Skin  return  versus  time 

8  -  Skin  return  histogram 

9  -  Azimuth  or  elevation  angle  error 
10  -  Psuedo-monopulse  error  signal. 

Not  all  plot  types  are  represented,  but  at  some  time  during  algorithm  development 
they  were  all  examined. 

Figures  5-1  and  5-2  give  baseline  results  for  COSRO  when  no  ECM  is  applied. 
Figures  5-3  through  5-9  are  sample  runs  for  the  parametric  controller  (see  Section 
2.2.1).  Figures  5-3  through  5-6  give  results  for  the  adaptive  controller. 

Figures  5-7  and  5-8  give  results  for  the  nonadaptive  controller  when  differing 
a  priori  information  is  given.  In  Figure  5-7  the  parameter  A,  which  has  a 
value  500  for  this  radar  model,  is  given  correctly  to  the  controller. 


L 


t 


i 


♦ 


i 


In  contrast,  Figure  5-8  shows  the  angle  error  when  the  parameter  A  is 
given  incorrectly  as  1250.  The  adaptive  controller,  assuming  no  prior 
knowledge  of  A  (set  equal  to  zero),  takes  longer  to  settle  to  near  the 
reference  angle,  as  expected;  see  Figure  5-9. 

Figures  5-10  through  5-20  are  sample  runs  for  the  input/output  controllers 
(see  Section  2.2.2).  Figures  5-10  through  5-12  are  samples  for  the  three 
different  controllers,  certainty  equivalence  (CE),  cautious,  and  dual,  all 
with  real-time  least  squares  (RTLS)  identification.  An  initial  identification- 
only  phase  is  done  on  all  runs.  Figures  5-13  through  5-18  show  comparisons 
of  differing  techniques  with  a  nonscintillating  target,  so  only  controller 
errors  are  shown.  Figure  5-13  is  an  open-loop  (step)  control  sample.  Figures 
5-14  through  5-16  are  the  three  controllers  with  RTLS,  while  Figures  5-17 
and  5-18  are  the  cautious  and  dual  controllers  with  recursive  instrumental 
variable  (RIV)  identification.  (The  CE  controller  was  virtually  identical  to 
the  cautious  controller.)  Figures  5-19  and  5-20  show  two  sample  runs  for 
the  controller  that  destabilizes  the  radar  at  a  controlled  rate.  Cautious 
control  with  RTLS  was  used. 


♦ 
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Table  5-1.  Sample  Runs  Against  COSRO 


FIGURE 

PLOT  TYPE 

REFERENCE  ANGLE 

CONTROLLER 

5-1 

7 

No  ECM 

5-2 

1 

- 

No  ECM 

5-3 

5 

0.3  BW 

Adaptive  parametric 

5-4 

1 

0.3  BW 

Adaptive  parametric 

5-5 

3 

0.3  BW’ 

Adaptive  parametric 

5-6 

6 

0.3  BW 

Adaptive  parametric 

5-7 

1 

0.3  BW 

Nonadaptive  parametric 
(see  Note  1) 

5-8 

1 

0.3  BW 

Nonadaptive  parametric 
(see  Note  2) 

5-9 

1 

0.3  BW 

Adaptive  parametric 
(see  Note  3) 

5-10 

1 

10  mrad 

CE  with  RTLS 
(see  Note  4) 

5-11 

1 

10  mrad 

Cautious  with  RTLS 
(see  Note  4) 

5-12 

1 

10  mrad 

Dual  with  RTLS 
(see  Note  4) 

5-13 

1 

- 

Open-loop  (see  Note  5) 

5-14 

1 

10  mrad 

CE  with  RTLS 
(see  Notes  4  and  5) 

5-15 

1 

10  mrad 

Cautious  with  RTLS 
(see  Notes  4  and  5) 

5-16 

1 

10  mrad 

Dual  with  RTLS 
(see  Notes  4  and  5) 

5-17 

1 

10  mrad 

Cautious  with  RIV 
(see  Notes  4  and  5) 

5-18 

1 

10  mrad 

Dual  with  RIV 
(see  Notes  4  and  5) 

5-19 

1 

Destabil ization 

Cautious  with  RTLS 
(see  Note  6) 

5-20 


1 


Destabilization 


Cautious  with  RTLS 
(see  Note  6) 


Motes 


1  - 

A  =  500, 

A=  500. 

2  - 

A  =  500, 

A  =  1250. 

3  - 

A-  500, 

o 

II 

<< 

4  -  Identification  only  until  pulse  number  100. 

5  -  Target  not  scintillating. 

6  -  Identification  only  until  pulse  number  200. 
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Figure  5-1 .  Skin  Return  -  No  ECM  -  COSRO 


5-3.  ECM  Signal  -  Adaptive  Parametric  Controller  -  COSRO 


Error  -  Adaptive  Parametric  Controller  -  COSRO 


Figure  5-5.  Angle  Error  -  Adpative  Parametric  Controller  -  COSRO 


Figure  5-6.  ECM  Histogram  -  Adpative  Parametric  Controller  -  COSRO 


Figure  5-8.  Angle  Error  -  Nonadaptive  Parametric  Controller,  A=  1250  -  COSRO 


Adaptive  Parametric  Controller,  A=  0  -  COSRO 


IS  IS  PL0T1S  NEXT  OHC.I2 


Figure  5-10.  Angle  Error  -  CE/RTLS  Controller  -  COSRO 


Figure  5-11.  AngLe  Error  -  Caut ious/RTLS  Controller  -  COSKO 
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Figure  5-12.  Angle  Krror  -  Dual/RTLS  Controller  -  COSRO 


g!*3H0  1X3M  SIIOI^  SI  S1M1 


Figure  5-13.  i\ngie  Error  -  Open-Loop  Controller  -  COSRO 


Figure  5-14.  Angle  Error  -  CE/KT1.S  Controller  -  COSRO 


Figure  5- II.  Angle  Error  -  Cant  Ions/ RTI.S  Gint.ro  I  lor  -  COSRO 


2i'3M0  J.X3N  SUOId  SI  SlHi. 


re  5-16.  Angle  Error  -  IXial/RTLS  Controller  -  CGSRO 


Figure  5-18.  Angle  Error  -  Dual/RTV  Controller  -  COSRO 


COSRO 


5-20.  Angle  Error  -  Cautious/RTLS  Controller  -  C0SR< 


5.2  SAMPLE  RUNS  FOR  PSUEDO-MONOPULSE 


Table  5-2  describes  the  sample  runs  done  against  the  psuedo -monopulse 
radar  system,  with  Figures  5-21  through  5-26  giving  outputs.  The  plot  types 
are  given  in  Sect  ion  5.1. 

Figure  5-21  is  the  baseline  result  for  psuedo -monopulse  when  no  ECM  is 
applied.  Figures  5-22  through  5-26  are  sample  runs  for  the  parametric 
controller  (see  Section  3.2.1);  no  runs  are  presented  for  the  input/output 
controller  for  this  radar  system.  The  angle  errors  for  the  azimuth  and 
elevation  is  given  separately  in  Figures  5-23  and  5-24,  respectively.  The 
90°  phase  shift  indicates  the  "orbiting"  behavior  of  the  pointing  angle. 

The  error  signal  at  the  output  of  the  summing  junction  (see  Figure  3-2)  is 
given  in  Figure  5-26  to  indicate  that  the  ECM  technique  does  indeed  introduce 
an  error  signal  that  is  not  subtracted  out. 


Table  5-2  Sample  Runs  Against  Psuedo-monopulse 


FIGURE 

PLOT  TYPE 

REFERENCE  ANGLE 

CONTROLLER 

5-21 

3 

0.3 

BW 

No  ECM 

5-22 

1 

0.3 

BW 

Adaptive  parametric 

5-23 

9 

(Azimuth) 

0.3 

BW 

Adaptive  parametric 

5-24 

9 

(Elevation) 

0.3 

BW 

Adaptive  parametric 

5-25 

5 

0.3 

BW 

Adaptive  parametric 

5-26 

10 

0.3 

BW 

Adaptive  parametric 

Error  -  No  ECM  -  Psuedo -Monopul  sc 


5-22.  Angle  Error  -  Adaptive  Parametric  Controller  -  Psuedo-Monopulse 


o 


Figure  'j-24 .  Elevation  Angle  Error  -  Adapt 


Figure  5-26.  Error  Signal  -  Adaptive  Parametric  Controller  -  Psuedo-Monopu  1  so 


5.3  SAMPLE  RUNS  FOR  MONOPULSE 


Table  5-3  describes  the  sample  runs  done  against  the  monopulse  radar 
system,  with  Figures  5-27  through  5-34  giving  outputs.  The  plot  types 
are  given  in  Section  5.1. 

Figures  5-27  through  5-34  are  sample  runs  for  the  input/output  controller 
(see  Section  4.2);  no  runs  are  presented  for  the  parametric  controller  for 
this  radar  system  as  the  controller  has  not  been  designed.  Figures  5-27 
through  5-30  are  using  the  cautious  controller  with  RTLS.  By  this  time 
it  was  decided  that  due  to  performance  and  ease  of  implementation,  this 
controller  is  probably  best  in  most  applications.  Angle  errors  and  ECM 
signals  are  given,  both  in  the  transition  region  and  in  steady-state.  A 
new  controller,  the  input  matching  controller  [1],  was  found  to  perform 
quite  well  for  this  system.  Its  outputs  are  given  in  Figures  5-31  through 
5-34.  This  controller  does  not  perform  identification  as  the  others  do, 
so  there  is  no  breakdown  of  time  into  these  two  phases. 


Table  5-3.  Sample  Runs  Against  Monpulse 


FIGURE 

PLOT  TYPE 

REFERENCE  ANGLE 

CONTROLLER 

5-2  7 

1 

0.3  BW 

Cautious  with  RTLS 
(see  Note  1) 

5-28 

5 

0.3  BW 

Cautious  with  RTLS 
(see  Note  1) 

5-29 

3 

0.3  BW 

Cautious  with  RTLS 
(see  Note  1) 

5-30 

5 

0.3  BW 

Cautious  with  RTLS 
(see  Note  1) 

5-31 

1 

0.3  BW 

Input  Matching 

5-32 

5 

0.3  BW 

Input  Matching 

5-33 

3 

0.3  BW 

Input  Matching 

5-34 

5 

0.3  BW 

Input  Matching 

Note 

1  -  Identification  only  until  pulse  number  200. 


5-32 


5-34 


Figure  5-28.  ECM  Signal  -  Caut ious/RTLS  Controller  -  Monopulse 


Figure  5-29.  Angle  Error  -  Cautious/ 


Figure  5-30.  ECM  Signal  -  Caut  ious/RTI.S  Controller  -  Monopulse 


Figure  5 


Figure  5-32.  ECM  Signal  -  Input  Matching  Controller  -  Monopulse 


Figure  5-33.  Angle  Error  -  Input  Matching  Controller  -  Monopulse 


Matching  Controller  -  Monopulse 
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